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Abstract 

The effects of crack inclination angle and crack length on the 
through-thickness mode I (KI), mode II (KII), and effective 
(. K e ff - y] Kf + Kfj) stress intensity factors (SIFs) have been 
analyzed by using three dimensional finite element analysis 
(3D FEA). Edge crack in semi circular bend specimen (SCB) 
was utilized in this investigation. The mode of mixity (Me) 
values are equal to 1 (pure mode I), 0.75, 0.5, 0.25, and 0 
(pure mode II). The crack length ratio, crack length/specimen 
radius (a/R), ranging from 0.1 to 0.7 by step equal toO.l has 
been studied. In SCB specimen, the mode I geometry 
correction factor (YI) decreased by increasing the crack 
length for all values of mode of mixity, while, in the case of 
Me = 1 & 0.75 YI reached its minimum value at a/R « 0.3 then 
YI increased by increasing the crack length. For all values of 
Me, Keff increased with the increasing crack length. 
However, the increment of increment in Keff for a/R < 0.3 is 
lower than that for a/R > 0.3. 
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Introduction 

Disc-type specimens are simple in geometry and have 
many advantages in terms of specimen preparation, 
testing and analysis. These test specimens have been 
used frequently to investigate mixed mode crack 
growth of rock materials, concrete, biomaterials, and 
other material [Fowell (1990, 1993), Lim (1993]. Semi- 
circular specimens, SCB, (Fig. 1) under three point 
bending were used for mixed Mode I-II fracture 
toughness calculation [Khan (2000), Lim (1994), 
Atkinson (1982)]. Depending on the crack length, its 
orientation with respect to the loading direction, and 
the distance between the supports, a variety of mixed- 
mode failure patterns can be achieved. For pure Mode- 
1, the crack is aligned parallel to the direction of the 
load and lies below the loading point. However, for 


pure Mode-II, a slight misalignment of the specimen 
and/or crack orientation may induce a component of 
Mode-I loading and pure Mode-II cannot be achieved 
[Khan (2000)]. The SCB specimen can give comparable 
results only when the proper span to diameter ratio is 
used. The specimen geometry requirement for a valid 
and representative fracture toughness value or stress 
intensity factor or other fracture properties of a rock 
material has been a matter of controversy among 
researchers. Crack length seems to be a more sensitive 
factor than specimen thickness on SIFs [Whittaker 
(1992)]. According to Lim [Lim (1994)], the SIFs 
become very sensitive at large a/R values for SCB 
specimens.The effect of specimen thickness and crack 
length on the variations of mode I and mode II stress 
intensity factors (SIFs) was analyzed either 
experimentally or using three dimensional finite 
element analysis (3D FEA) in Center cracked circular 
disc specimen (CCCD) and SCB specimen by Sallam et 
al. [Sallamand Abd-Elhady (2012), Mubarakiet al. 
(2013)]. They found that, the normalized mode I SIF 
increased at the mid plane of specimen by increasing 
the thickness of the specimen. However,at the 
specimen surface, this value decreased by increasing 
the specimen thickness. 

The present study investigated the effects of the 
sample geometry in terms of crack length and crack 
inclination angle on mode I and mode II stress 
intensity factors (SIFs) of edge crack in semi circular 
bend specimen (SCB), with the aim to make 
improvements in the testing techniques for fracture 
testing of brittle materials. 

Numerical Analysis 

The basic dimensions of SCB specimens R (specimen 
radius), and B (specimen thickness), were equal to 75 
mm and 22.5 respectively. SCB specimen is placed on 
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two bottom supports of distance 2S and the ratio of 
S/R was 0.43. The crack length ratio, a/R, ranged from 
0.1 to 0.7 by step equal toO.l. Aliha [Aliha 
(2010)]suggested five values for mode mixity by 
parameter called mixity parameter, M e , 

2 /K,\ 

M e — — arctanl— — ) (1) 

7T Vtt/// 

The values of M e varied through 1 (pure mode I), 0.75, 
0.5, 0.25and 0 for pure mode II. Ayatollahi and Aliha 
[Ayatollahi (2007)] conducted extensive finite element 
analyses and showedthat,this mixity parameter was 
valid by change the crack inclination angle, as 0, 
18.5, 33, 42.5 and 50°. 



SPECIMENS SUBJECTED TO MIXED MODE I/II LOADING 


In the present analysis, the mode I and mode II 
normalized stress intensity factors are denoted as 
Y/and Yu, respectively, and it can be deduced from 
Refs. [Hutar (2010), Lim (1993)] that the general 
formula for normalized stress intensity factor Yi, is 
defined as: 


Where: 


4 RtKi 
y. = - 

1 Pyfcur 


i = /, II 


( 2 ) 


Ki = mode I stress intensity factor 
Kn = mode II stress intensity factor 
P = applied load 

a, R, and t = crack length, radius of specimen, and half 
specimen thickness = £>/ 2, respectively. 

the mechanical properties of the specimen were as 
follows: modulus of elasticity, E = 54 Gpa, and 
poisson's ratio, v = 0.276. The specimen material is 
homogeneous, isotropic and elastic, the applied load, 
P was equal 50 kN as shown in fig. 1. the plane x- 
y(plane z = 0) is the mid plane of specimen and two 
specimen surface are z = t and z = -f, respectively, in 


order to calculate stress intensity factors of the 
samples with different geometries, numerical 
computations were carried out. the package programs 
used in this study were ABAQUS. ABAQUS program, 
a finite element program, was used in the present 
work and all the numerical studies were conducted 
with3D models. Disc type specimens require 3D 
modeling for stress analysis and stress intensity factor 
computations. For 3D,modeling a choice had to be 
done among ABAQUS. Around 20 planar layers are 
divided through the thickness of the SCB specimen. 
Within each layer, the size of element, C3D8 "8-node 
linear brick element", decreases gradually with 
distance from the crack tip decreasing. That is mean 
FE meshes in the neighborhood of the crack tip are 
much denser. Figure 2 shows a typical example of the 
present idealization. 



FIG. 2.TYPICAL FINITE ELEMENT MESH IDEALIZATION OF 
THE SCB SPECIMEN. 


To determine if this program produces accurate results 
for stress intensity factor calculations comparable to 
the result of other references [Ayatollahi (2007)], 
simple models with known analytical results for 
model verification were first employed. After 
verification studies, the program was determined to be 
appropriate for further stress intensity factor 
computations. In the present work, the domain 
integral methodiscommonly used to extract stress 
intensity factors (SIFs) [Nakamura (1989, 1991), Gosz 
(1998, 2002)]. In a finite element, model SIF can be 
thought of as the virtual motion of a block of material 
surrounding each node along thecrack line. Each such 
block is defined by contours: each contour is a ring of 
elements completely surrounding the nodes along the 
crack line from one crack face to the opposite crack 
face. These rings of elements are defined recursively to 
surround all previous contours. ABAQUS/ Standard 
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automatically finds the elements that form each ring 
from the regions given as the crack-line definition. 
Each contour provides an evaluation of the contour 
integral. Using the divergence theorem, the contour 
integral can be expanded into a volume integral, over 
a finite domain surrounding the crack. This domain 
integral method is used to evaluate contour integrals 
in ABAQUS/ Standard. The method is quite robust in 
the sense that accurate contour integral estimates are 
usually obtained even with quite coarse meshes; 
because the integral is taken over a domain of 
elements surrounding the crack, errors in local 
solution parameters have less effect on the evaluated 
quantities such as the stress intensity factors. The 
stress intensity factors KI, KII and Kill (mode I, mode 
II and mode III SIF) are usually used in linear elastic 
fracture mechanics to characterize the local crack- 
tip/crack-line stress and displacement fields. They are 
related to the energy release rate (the J- integral). The 
energy release rate is calculated directly in ABAQUS/ 
Standard. 

Results 

Figure 3 shows the results obtained from finite 
element analysis in the present work for a/R = 0.3 and 
B/R = 0.1 and S/R = 0.43 and the results given by 
Ayatollahi andAliha [Ayatollahi (2007)] at the surface 
of the SCB specimen, z = t. There is a good agreement 
between the two sets of results and it can be 
considered as a validation for the present analysis. 


value of z/f,meaning that the value of Yi at surface 
point of SCB specimen was not the highest value on 
the crack front. By increasing the crack length, the 
value of Yi decreases to reach a minimum value then it 
increases that for ATequal to 1 and 0.75 only that 
agrees with Lim [Lim (1994)]. Otherwise, for AT = 0.5 
and 0.25, the value of Yi decreases by increasing the 
value of crack length. 



0 0.1 0.2 0.3 0.4 0.5 

a/R 


0.1 0.2 0.3 0.4 

a/R 


FIG. 4 THE EFFECT OF SCB SPECIMEN CRACK LENGTH ON 
THE NORMALIZED MODE I SIF FOR DIFFERENT THE SITE ON 
THE CRACK FRONT z/t AND DIFFERENT MODE MIXITY A/I' (A) 
AT = 1, (B) M' = 0.75, (C) AT = 0.5 AND (D) AT = 0.25. 
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FIG. 3. NORMALIZED MODE I AND MODE SIFS CALCULATED 
IN PRESENT WORK COMPARED WITH THOSE PRESENTED BY 
AYATOLLAHI AND ALIHA [Ayatollahi (2007)]FOR THE SCB 
SPECIMEN ( a/R = 0.3, B/R = 0.1 AND S/R = 0.43) 



Figure 4 depicts the relation between normalized 
mode I SIF and crack length at different site at crack 
front, z/t, for SCB specimen with different mode 
mixity. The value of Yi increases by decreasing the 


FIG. 5 THE EFFECT OF SCB SPECIMEN CRACK LENGTH ON 
THE NORMALIZED MODE II SIF FOR DIFFERENT THE SITE ON 
THE CRACK FRONT z/t AND DIFFERENT MODE MIXITY AT (A) 
At' = 1, (B) At" = 0.75, (C) AT' = 0.5 AND (D) At' = 0.25. 
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A little effect of z/f on the value of normalized mode II 
SIF, Yn, with different crack length can be shown in 
Fig. 5. The value of Yu decreases by increasing the 
value of crack length to reach the minimum value then 
increases by increasing the crack length that for AT = 
0.5, 0.25 and 0. However, for AT = o.75 the value of Yn 
increases by increasing the crack length. 


Conclusions 

The present numerical analysis revealed the following 
conclusions: 

1- For a/R< 0.3, the mode I geometry correction 
factor ( Yi ) decreased by increasing the crack 
length for all values of mode of mixity. 


Figure 6 shows the effect of crack length, a/R on the 
effective stress intensity factor. Keg, {K ef f = yjKf + K? t ) 
through the SCB specimen crack front. The change of 
the effective stress intensity factor, Keff, ( K eff - 
V Kf + Kf t ) through the SCB specimen crack front can 
be neglected. By increasing the crack length, the values 
of fu# increase. In the case of pure mode I and II, the 
effective stress intensity factor rises by increasing the 
crack length. 




a/R 

FIG. 6 THE EFFECT OF SCB SPECIMEN CRACK LENGTH ON 
THE SIF EFFECTIVE FOR DIFFERENT THE SITE ON THE CRACK 
FRONT z/f AND DIFFERENT MODE MIXITY M c (A) M‘ = 1, (B) M‘ 
= 0.75, (C) M‘ = 0.5 AND (D) M‘ = 0.25. 


2- For a/R> 0.3, Y/increased by increasing the crack 
length for AT = 1 and 0.75, while, Yi continued to 
decrease with the increasing crack length for 
Me= 0.5and 0.25. 

3- For a/R> 0.3, Y/zincreased by increasing the crack 
length for all values of Me. 

4- For all values of /IT, Keff increased with the 
increasing crack length. However, the increment 
of increment in Keff for a/R <0.3 was lower than 
that for a/R > 0.3. 
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